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Study of jets produced in association with a W boson in pp collisions 
at ^/s = 7 TeV with the ATLAS detector 



(N 

o 

(N 

>> 
c3 



X 



> 

(N 



O 



X 



The ATLAS Collaboration* 
(Dated: January 24, 2012) 

We report a study of final states containing a W boson and hadronic jets, produced in proton- 
proton collisions at a center-of-mass energy of 7 TeV. The data were collected with the ATLAS 
detector at the CERN LHC and comprise the full 2010 data sample of 36 pb _1 . Cross sections are 
determined using both the electron and muon decay modes of the W boson and are presented as 
a function of inclusive jet multiplicity, iVj et , for up to five jets. At each multiplicity, cross sections 
are presented as a function of jet transverse momentum, the scalar sum of the transverse momenta 
of the charged lepton, missing transverse momentum, and all jets, the invariant mass spectra of 
jets, and the rapidity distributions of various combinations of leptons and final-state jets. The 
results, corrected for all detector effects and for all backgrounds such as diboson and top quark pair 
production, are compared with particle-level predictions from perturbative QCD. Leading-order 
multiparton event generators, normalized to the NNLO total cross section for inclusive W-boson 
production, describe the data reasonably well for all measured inclusive jet multiplicities. Next-to- 
leading-order calculations from MCFM, studied here for iVjet < 2, and BlackHat-Sherpa, studied 
here for JVjet < 4, are found to be mostly in good agreement with the data. 

PACS numbers: 12.38.Qk, 13.85.Hd, 13.85.Qk, 13.87.Ce 



I. INTRODUCTION 

The study of massive vector boson production in as- 
sociation with one or more jets is an important test 
of quantum chromodynamics (QCD). These final states 
are also a significant background to studies of Standard 
Model processes such as ti, diboson, and single-top pro- 
duction, as well as to searches for the Higgs boson and 
for physics beyond the Standard Model. Thus, measure- 
ments of the cross section and kinematic properties, and 
comparisons with theoretical predictions, arc of signifi- 
cant interest. Measurements of VL+jets production in 
proton-antiproton collisions at s/s = 1.96 TeV have been 
reported by the CDF and DO Collaborations [1, 2] and 
for y/s = 7 TeV proton-proton collisions by the CMS 
Collaboration [3]. Measurements of jets produced in as- 
sociation with a Z boson were also performed using pp 
collisions at y/s = 1.96 TeV [4-6] and pp collisions at 
yfs = 7 TeV [3, 7]. The study presented here is comple- 
mentary to the measurement of the transverse momen- 
tum distribution of W bosons conducted by the ATLAS 
Collaboration [8]. 

This paper reports a measurement at the CERN Large 
Hadron Collider (LHC) of the VL+jets cross section for 
proton-proton (pp) collisions at a center-of-mass energy 
(y/s) of 7 TeV, using the ATLAS detector. The measure- 
ment is based on the full 2010 data sample, corresponding 
to an integrated luminosity of approximately 36 pb _1 . 
It is an extension of an earlier ATLAS measurement of 
both the electron and muon decay modes of the W boson 
based on 1.3 pb _1 [9]. Compared to the earlier result, 
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uncertainties in both the jet energy scale and luminos- 
ity are reduced, acceptance for the jets is expanded, and 
event reconstruction and simulation are improved. The 
improved reconstruction brings better alignment of the 
detector systems and reduction of backgrounds in the 
electron channel. 

The results have been corrected for all known detec- 
tor effects and are quoted in a specific range of jet and 
lepton kinematics, fully covered by the detector accep- 
tance. This avoids modcl-dcpcndcnt extrapolations and 
facilitates comparisons with theoretical predictions. The- 
oretical calculations at ncxt-to- leading order (NLO) in 
perturbative QCD (pQCD) have been computed inclu- 
sively for up to four jets [10, 11] and are compared with 
the data. 



II. THE ATLAS DETECTOR 

ATLAS uses a right-handed coordinate system with 
its origin at the nominal pp interaction point (IP) in the 
center of the detector and the z-axis along the beam 
pipe. The x-axis points from the IP to the center of 
the LHC ring, and the y-axis points upward. Cylindri- 
cal coordinates (r, <fi) are used in the transverse plane, <j> 
being the azimuthal angle around the beam pipe. The 
pseudorapidity is defined in terms of the polar angle 
6 as 77 = — lntan(#/2) and the rapidity is defined as 
y = \n[(E + p z )/(E — p z )}/2. The separation between 
final state particles is defined as AR = \J (Ay) 2 + (A0) 2 
and is Lorentz invariant under boosts along the z-axis. 

The ATLAS detector [12, 13] consists of an inner 
tracking system (inner detector, or ID) surrounded by 
a thin superconducting solenoid providing a 2T mag- 
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netic field, electromagnetic and hadronic calorimeters, 
and a muon spectrometer (MS). The ID consists of pixel 
and silicon microstrip detectors, surrounded by a tran- 
sition radiation tracker. The electromagnetic calorime- 
ter is a liquid-argon and lead detector, split into bar- 
rel (\t)\ < 1.475) and endcap (1.375 < \r]\ < 3.2) re- 
gions. Hadron calorimctry is based on two different 
detector technologies. The barrel (\r/\ < 0.8) and ex- 
tended barrel (0.8 < < 1.7) calorimeters arc com- 
posed of scintillator and steel, while the hadronic end- 
cap calorimeters (1.5 < |7y| < 3.2) utilize liquid-argon 
and copper. The forward calorimeters (3.1 < < 4.9) 
are instrumented with liquid-argon/copper and liquid- 
argon/tungsten, providing electromagnetic and hadronic 
energy measurements, respectively. The MS is based 
on three large superconducting toroids arranged with an 
eight-fold azimuthal coil symmetry around the calorime- 
ters, and a system of three stations of chambers for trig- 
gering and for precise track measurements. 

III. DATA AND ONLINE EVENT SELECTION 

The data for this analysis were collected during LHC 
operation in 2010 with proton-proton interactions at a 
center-of-mass energy of 7 TeV. The collisions occurred 
within pairs of bunches of up to ~ 1.1 x 10 11 protons 
per bunch. The bunches were configured in trains with a 
time separation between bunches of 150 ns and a longer 
separation between trains. Data were collected with up 
to 348 colliding bunch pairs per beam revolution. This 
configuration led to a peak instantaneous luminosity of 
up to 2.1 x 10 32 cm _2 s _1 that corresponds to an aver- 
age of 3.8 inelastic collisions per bunch crossing. Typical 
values were lower as the luminosity degraded during the 
data taking fills which lasted up to 20 hours. On aver- 
age, the data contain 2.1 inelastic collisions per bunch 
crossing. 

Application of beam, detector, and data-quality re- 
quirements resulted in a total integrated luminosity of 
36 pb _1 . The uncertainty on the luminosity is 3.4% [14, 
15]. The integrated luminosities for the data samples as- 
sociated with the electron and muon decay modes of the 
W boson were calculated separately and differ by 1.7%. 

Events were selected online if they satisfied either the 
electron or muon criteria described below. Criteria for 
electron and muon identification, as well as for event se- 
lection, followed closely those of the previous 1.3 pb _1 
W+jets cross section analysis [9]. 

For this analysis, the following kinematic requirements 
were imposed on events in order to enter the selected 
sample: 

• p T > 20 GcV (£ = electron or muon), 

• \ri e \ < 2.47 (except 1.37 < \if\ < 1.52) or \?f\ <2.4, 

• _E™ 18S > 25 GcV (missing transverse momentum), 

• m T (W) > 40 GcV, 



• ^ > 30 GeV > 

• \yi ct \ < 4.4 and AR(£,jet) > 0.5. 

These selection criteria differ slightly from the fiducial ac- 
ceptance to which measured cross sections are finally cor- 
rected, which is described in Section VF. The transverse 
momenta of the leptons and neutrinos from W — > ev and 
W — > [iv decays are denoted as p T and p T , respectively. 
The transverse momentum of the neutrino is determined 
as E?f- lss , the missing transverse momentum, from the 
requirement that the total transverse momentum of all 
final-state particles is a zero vector. The calculation of 
E™ lss and the transverse mass of the W, m^(W), are 
discussed later in Section VB. 

All measured cross sections are corrected for any de- 
tection losses within these regions. The lower bound 
p^ > 30 GeV is chosen to facilitate comparisons with 
other experiments and with next-to-leading-ordcr QCD 
predictions. Appendix A shows analogous results with 
Pr* > 20 GeV in order to facilitate validation of the QCD 
description in Monte Carlo generators and future theo- 
retical developments in this area. 



A. Electron selection 

In the electron channel, events were selected online us- 
ing two different triggers depending on the instantaneous 
luminosity. The tighter trigger requirement corresponds 
to 99.1% of the data and is a subset of the looser one. 
It required the presence of at least one electromagnetic 
cluster in the calorimeter with transverse energy above 
15 GeV in the region of |?7| < 2.5. The final selection 
requirements were applied by the online event filter [12] 
and the kinematic variables correspond closely to those 
in the offline analysis described in Section VC. 

The impact of the trigger efficiency was small for elec- 
trons with Et > 20 GcV, as required in this analysis. The 
efficiency was measured using Z — > ee decays identified 
in the experimental data. It was found to be 99.0 ± 0.5% 
and constant over the full kinematic region of this mea- 
surement [16, 17]. 



B. Muon selection 

In the muon channel, events were selected online using 
a trigger that required the presence of a muon candidate 
reconstructed in both the muon spectrometer and inner 
detector, consistent with having originated from the in- 
teraction region. The candidate was required to have 
Pt > 10 GeV or px > 13 GeV (depending on the data- 
taking period) and |f?| < 2.4. The higher threshold was 
used to collect most of the data. As in the electron case, 
these requirements were imposed in the online event filter 
and were less stringent than those applied offline. The 
offline selection is documented later in Section VD. The 
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average trigger efficiency was measured to be ^85% in- 
cluding the reduced geometrical acceptance in the central 
region. 



IV. SIMULATED EVENT SAMPLES 

Simulated event samples were used for most back- 
ground estimates, for the correction of the signal yield 
for detector effects and for comparisons of results to the- 
oretical expectations. The detector simulation [31] was 
performed using GEANT4 [32]. The simulated event 
samples are summarized in Table I for signal simulations 
and Table II for the background simulations. The Alp- 
GEN and Mc@Nlo samples were interfaced to Herwig 
for parton shower and fragmentation processes and to 
Jimmy v4.31 [33] for underlying event simulation. Sim- 
ilarly, Jimmy was used for the underlying event simula- 
tion in the diboson samples produced with Herwig. The 
AcerMC ti samples were showered with Pythia where 
the default settings for initial state radiation (ISR) and 
final state radiation (FSR) were altered [34] . The param- 
eterization of the factorization scale used for the matrix 
element (ME) calculation in the Alpgen samples was 

chosen to be Q\ = m v + ^ (Pt)s where my is the 

partons 

mass of a W or Z boson and the decay products of the 
boson are not included in the sum [21]. The parton-jet 
matching was performed at = 20 GeV with the MLM 
matching scheme [35] using jets from the cone clustering 
algorithm with R = 0.7. The default renormalization and 
factorization scales were used in the Sherpa samples and 
the parton-jet matching was performed at = 30 GeV 
using the CKKW matching scheme [36, 37]. Parton den- 
sity functions (PDFs) were: CTEQ6L1 [38] for the Alp- 
gen samples and the parton showering and underlying 
event in the PowHeg samples interfaced to Pythia; 
MRST 2007 LO* [39] for Pythia, AcerMC, and the 
diboson samples; and CTEQ6.6M [40] for Mc@Nlo, 
Sherpa, and the NLO matrix element calculations in 
PowHeg. The radiation of photons from charged lcp- 
tons was treated in Herwig and Pythia using Photos 
v2.15.4 [41]. Tauola vl.0.2 [42] was used for r lepton de- 
cays. The underlying event tunes were the ATLAS MC10 
tunes: ATLAS underlying event tune #1 (AUET1) [43] 
for the Herwig, Alpgen, and Mc@Nlo samples; AT- 
LAS minimum bias 1 (AMBT1) [44] for Pythia, Ac- 
erMC, and PowHeg samples. These two tunes were 
derived using pp collisions at ^Js = 7 TeV produced at 
the LHC. The samples generated with Sherpa used the 
default underlying event tune determined from lower en- 
ergy measurements and pp data from the LHC. 

Samples were generated with minimum bias interac- 
tions overlaid on the hard-scattering event to account for 
the multiple pp interactions in the same beam crossing 
(pile- up) . The minimum bias interactions were simulated 
with Pythia with the AMBT1 tune. These samples were 



then re- weighted so the distribution of the number of pri- 
mary vertices matched that of the data. 



V. OFFLINE EVENT ANALYSIS 

Events were selected if they satisfied the criteria de- 
scribed above and had at least one interaction vertex with 
three or more associated charged particle tracks, located 
within 200 mm in z from the center of the detector. For 
these data the luminous region had a typical RMS size of 
~ 60 mm in z. The position resolution of reconstructed 
vertices along z was ~ 0.1 mm for a vertex with 10 tracks. 
For the sample of events passing the single-lepton trig- 
ger the mean number of interaction vertices was 2.1 per 
event. The primary vertex was taken as the one with the 
largest £p T of associated tracks. Events with significant 
noise in the calorimeters, cosmic rays, and beam-induced 
background were rejected [45]. 

A. Jet selection 

Jets were reconstructed from energy observed in the 
calorimeter cells using the anti-fc t algorithm [46] with a 
radius parameter R = 0.4 [45]. Since the volume of in- 
dividual cells is small compared to the volume of the 
electromagnetic and hadronic energy showers, cells were 
grouped into clusters depending on their signal size rela- 
tive to noise [47] . These clusters formed the input to the 
jet reconstruction. Since a jet involves many clusters a 
mass can be calculated and the jet rapidity rather than 
pseudorapidity was determined. 

To account for the difference in calorimeter response 
between electrons and hadrons of the same energy, and 
to correct for other experimental effects, a pt and 77- 
dependent factor, derived from simulated events, was 
applied to each jet to provide an average energy scale 
correction [45]. Jets were required to have a rapidity 
\y\ < 4.4 and pt > 30 GeV. To ensure a reliable energy 
measurement all jets within AR < 0.5 of an electron 
or muon (that passed the lepton identification require- 
ments) were explicitly not considered, regardless of the 
jet pt or rapidity, but the event itself was retained. Jets 
consistent with detector noise, cosmic rays, or beam halo 
were rejected [45]. The jet rejection requirement was 
more stringent than that applied to events. 

To suppress jets arising from additional pp interactions 
a parameter called the jet-vertex fraction (JVF) was cal- 
culated for each jet in the event. After associating tracks 
to jets by requiring AR < 0.4 between tracks and a jet, 
the JVF was computed for each jet as the scalar sum 
of pt of all associated tracks from the primary vertex 
divided by the total pr associated with that jet from 
all vertices. The JVF could not be calculated for jets 
which fell outside the fiducial tracking region (\r)\ < 2.5) 
or which had no matching tracks so these were assigned 
a value of —1 for accounting purposes. Only jets with 
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TABLE I. Samples of simulated signal events used in this analysis. The W samples are normalized to the inclusive NNLO 
cross section of 10.46 nb calculated with FEWZ [18] using the MSTW2008 PDF set [19]. For PYTHIA, the inclusive W sample 
is based on a 2 —¥ 1 matrix element merged with a 2 — > 2 matrix element and a leading-logarithmic parton shower. Details of 
PDF sets, final-state photon radiation, and underlying event tunes are given in the text. 



Physics process 


Generator 


W inclusive (W 


-+£v;£ = 


e,fi,r) 




Pythia 6.4.21 


[20] 


W + jets (W -> 


iu;£ = e 


fx; < JVparton 


<5) 


Alpgen 2.13 


[21] 


W + jets (V -> 


£u;£ = e 


Mi < JVparton 


<5) 


Sherpa 1.3.1 


[22] 



TABLE II. Samples of simulated background events used in this analysis. The Z+jets samples were normalized using the 
inclusive cross sections from FEWZ [18] code that utilized MSTW2008 PDF set [19]. The tt cross section is given at next-to- 
leading order (plus next-to-next-to-leading-log, NNLL). The dijet cross sections are given at leading order (LO) in pQCD. For 
these samples, the variable pt is the average px of the two outgoing partons from the hard-scattering process before modification 
by initial- and final-state radiation and the underlying event. Details of PDF sets, final-state photon radiation, and underlying 
event tunes are given in the text. 
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Alpgen 2.13 


[21] 
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NNLO 


[18] 


Z -^tt (m« > 60 GeV) 


Pythia 6.4.21 


[20] 


0.989 


NNLO 


[18] 


tt 


PowHeg-HVQ 












vl.01 patch 4 


[23] 


0.165 


NLO+NNLL 


[24] 


tt 


AcerMC 3.7 


[25] 


0.165 


NLO+NNLL 


[24] 


Single-top t — > Ivq (s-channel) 


McONlo 3.3.1 


[26, 27] 


4.3xl0~ 4 


NLO 


[28] 


Single-top t — > Ivq (t-channel) 


McONlo 3.3.1 


[26, 27] 


6.34xl0~ 3 


NLO 


[28] 


Single-top (Wt) 


McONlo 3.3.1 


[26, 29] 


13.1xl0~ :i 


NLO 


[28] 


WW 


Herwig 6.510 


[30] 


44.9xl0~ 3 


NLO 


[28] 


WZ (m z > 60 GeV) 


Herwig 6.510 


[30] 


18.5xl0~ 3 


NLO 


[28] 


ZZ (m z > 60 GeV) 


Herwig 6.510 


[30] 


5.96 xl0~ 3 


NLO 


[28] 


Dijet (/J channel, p T > 8 GeV, p T > 8 GeV) 


Pythia 6.4.21 


[20] 


10.6 xlO 6 


LO 


[20] 



the absolute value of the JVF smaller than 0.75 were re- 
jected so that jets with a JVF of —1 were kept. Figure 1 
shows the distribution of this parameter for all jets in the 
W — >• ev data and Monte Carlo event samples. The re- 
quirement on the JVF is most important for low px jets 
and for the data with high instantaneous luminosity. 

The pile-up collisions also add a uniform background of 
particles to the events and slightly increase the measured 
jet energies. The jet energy calibration factor described 
above contains a correction for this effect. 

No minimum separation AR was required between fi- 
nal state jets, but the measured jet response changed for 
separations less than AR < 0.5. This distortion in the 
response was corrected by the event reconstruction effi- 
ciency calculation and residual effects enter the estimated 
systematic uncertainties. 

After the application of all jet requirements, the effi- 
ciency for reconstructing jets was determined from sim- 
ulation to be ~ 97% for jets with px = 30 GeV, rising 
to close to 100% for jets above 80 GeV. The uncertain- 
ties in the jet energy scale and jet energy resolution were 
determined in separate studies [45] . The uncertainties in 
the jet energy scale were 2.5-14%, and depended on the 
77 and px of the jet. The uncertainty on the jet energy 
resolution was ~ 10% for each jet, relative to the nominal 
resolution which also varied with r\ and px- 



B. Missing transverse momentum and m,T(W) 

The calculation of missing transverse momentum 
(i5™ 18S ) and transverse mass of W bosons (i71t(W)) fol- 
lowed the prescription in Rcfs. [48] and [16]. i71t(W) was 
defined by the lepton and neutrino px and direction as 

m T (W) = ^/2p T p T (l - cos((/) e - (j> v )), where the (x,y) 

components of the neutrino momentum were taken to be 
the same as the corresponding i?™ lss components. E™ ss 
was calculated from the energy deposits in calorimeter 
cells inside three-dimensional clusters [47] . These clusters 
were then corrected to account for the different response 
to hadrons compared to electrons or photons, as well as 
dead material and out-of-cluster energy losses [49] . Only 
clusters within \rj\ < 4.5 were used. In the muon chan- 
nel, £^™ lss was corrected for the muon momentum and its 
energy deposit in the calorimeters. Events were required 
to have £ T liss > 25 GeV and m T (W) > 40 GeV. 



C. W —¥ ev + jets final state 

Electrons were required to pass the standard "tight" 
electron selection criteria [16, 17] with £x > 20 GeV and 
|ry| < 2.47. Electrons in the transition region between the 
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FIG. 1. Jet-vertex fraction (JVF) distribution for all jets in 
the W — > ev sample. The events at —1 correspond to jets 
where the JVF could not be calculated, while the peak near 
corresponds to jets from a secondary vertex. For the data 
99.1% of the jets pass the requirement that the absolute value 
of the JVF be greater than 0.75, while for the Monte Carlo 
sample this rate is 98.8%. 

barrel and endcap calorimeter (1.37 < \r\\ < 1.52) were 
rejected. 

To suppress multijet events containing non-isolated 
electrons such as those from semileptonic decays of 
hadrons containing charm and bottom quarks, a 
calorimeter-based isolation requirement was applied. 
The transverse energy within a cone of radius R — 0.2 
around the electron, corrected for contributions from the 
electron, was required to be less than 4 GcV. This isola- 
tion requirement is more than 96% efficient over all jet 
multiplicities for prompt electrons originating from de- 
cays of W bosons and reduces the non-isolated electron 
background by a factor of two. 

To remove backgrounds from Z — > ee decays, events 
were also rejected if there was a second electron passing 
the "medium" electron selection criteria [16, 17] and the 
same kinematic selections and isolation requirements as 
above. 



1. Electron channel background estimates 

The principal backgrounds in the electron channel arise 
from multijet QCD events, other leptonic decays of gauge 
bosons and, at higher jet multiplicities, ti production. 
The background from gauge bosons includes W — > rv 
where the r lepton decays to an electron and Z — > ee 
where one electron is not identified and hadronic en- 
ergy in the event is mismcasurcd. Leptonic ti decays 



(ti —> bbqq'ev), single-top events and diboson (WW, 
WZ, ZZ) processes were also evaluated. The number of 
leptonic background events surviving the above selection 
requirements was estimated with simulated event sam- 
ples that were introduced earlier in Section IV. Specifi- 
cally, Pythia was used for W —> tv and Z — > tt and 
Alpgen for the other vector boson samples. The sim- 
ulated leptonic background samples were normalized to 
the integrated luminosity of the data using the predicted 
cross sections shown in Table II. The ti background is 
discussed in more detail later in Section VE. 

The multijet background in the electron channel has 
two components, one where a light flavor jet passes the 
electron selection and additional energy mismeasurement 
results in large E!f lss , and the other where a bottom or 
charm hadron decays to an electron. The number of 
multijet background events was estimated by fitting, for 
each exclusive jet multiplicity, the E™ lss distribution in 
the data (without the E™ lss selection requirement) to a 
sum of two templates: one for the multijet background 
and another which included signal and the leptonic back- 
grounds. The fits determined the relative normalizations 
of the two templates for each exclusive jet multiplicity. 
The shapes for the second template were obtained from 
simulation and their relative normalization was fixed to 
the ratio of their predicted cross sections. 

The template for the multijet background was obtained 
from the data because the mechanisms by which a jet 
fakes an electron are difficult to simulate reliably. The 
template was derived by loosening some of the electron 
identification requirements. Two approaches were taken 
so their results could be compared. 

In the first, the requirements on shower shape in the 
calorimeter were relaxed. The "loose" electron identifi- 
cation criteria of Rcf. [16, 17] were applied to the shower 
shapes. The track-cluster matching requirements applied 
in the standard "tight" electron selection were still ap- 
plied but the remaining "tight" requirements with re- 
spect to the "medium" requirements were required to 
fail [16, 17]; the selection favors electron candidates from 
conversions or from charged hadrons overlapping electro- 
magnetic showers. 

In the second method, the requirement that a track 
matched the energy deposition in the calorimeter was 
relaxed and loose photon identification requirements were 
used instead of those of an electron. 

To suppress any residual signal contribution, the iso- 
lation requirement was also reversed in both methods. 
A large simulated dijet sample was used to verify that 
these requirements do not bias the E™ lss shape of the 
background templates. 

The results of the two methods were compared for each 
jet multiplicity and agreed within their statistical uncer- 
tainties. For the zero-jet bin they agreed to better than 
17% with respect to the total number of candidate back- 
ground events. Residual differences are included in the 
estimates of systematic uncertainty described below. The 
range of _B™ 1SS used to fit the templates was also varied to 
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FIG. 2. Result of the _B™ 1SS template fits used to obtain an estimate of the multijet background for W — > events, in bins 
of exclusive jet multiplicity. The data is shown with the statistical uncertainties only. In this case the multijet template 
was obtained with relaxed shower shape requirements, as described in the text. The data with > 5 jets are not used for 
measurements because of low event multiplicity and a poor signal-to-background ratio. The event multiplicity and the ratio 
were better for p^' > 20 GeV. W candidate events were required to have E™ lss > 25 GeV. 



estimate systematic effects. The first method was used to 
calculate the central values of the multijet backgrounds 
for the various jet multiplicities. 

The comparisons of the template fits to the E!f lss dis- 
tributions are shown in Fig. 2 for the first type of multijet 
template. Figure 3 shows the final tot(W / ) distributions 
in the various bins of inclusive jet multiplicity. 



2. Electron channel systematic uncertainties 

The systematic uncertainties for the electron channel 
are summarized in Table III. The calculation of uncer- 
tainty on the number of multijet background events was 
introduced in Section VC 1. 

The electron trigger efficiency was measured using 
Z — > ee events triggered by an object other than the 
electron under study (tag-and-probc method). A scale 
factor of 99.5 ± 0.5% relative to the value predicted by 
the Monte Carlo simulation was determined. The same 
event samples were used to determine the electron re- 
construction and identification efficiencies relative to the 
Monte Carlo prediction. The reconstruction efficiencies 



were consistent with the Monte Carlo values within a 
systematic uncertainty of 1.5%. Data-driven corrections 
to the simulated identification efficiencies were charac- 
terized by a two-dimensional matrix in r\ and £t- The 
Z — > ee events were also used to test the electron identi- 
fication efficiency for any dependence on accompanying 
jet activity and none was found. 

The measured electron energy scale and resolution 
were also studied with Z — > ee events. In the data, elec- 
tron energies were adjusted with an Ty-dcpcndcnt correc- 
tion with typical values of about 2% [17]. The electron 
energy resolution was similarly tested and adjusted in 
simulated events. The residual systematic uncertainties 
arc shown in Table III. 



D. W — > [iv + jets final state 

The muons were required to be reconstructed in both 
the ID and MS subsystems and to have pt > 20 GeV 
and I77I < 2.4. The ID track requirements were those 
of Ref. [16]. An ID-based muon isolation was applied 
which required a relative isolation of SpJjP/p^, < 0.1, us- 
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FIG. 3. Transverse mass distributions mT(W) for selected W — > ev events in bins of inclusive jet multiplicity. MC predictions 
for the signal and leptonic backgrounds are normalized to luminosity using (N)NLO cross sections and the multijet background 
is estimated from data (method I). 



TABLE III. Summary of relative systematic uncertainties as- 
sociated with the electron channel. 



Quantity 


Uncertainty 


Trigger efficiency 
Electron reconstruction 
Electron identification 
Electron energy scale 
Electron energy resolution 
Multijet QCD background 


~ 0.5% 
~ 1.5% 
2 - 8% a 
0.3- 1.6% a 
< 0.6% of the energy 
17 - 100% b ; difference between the 
two methods, see Section V C 1 



a V — Pt dependent 

b increased with jet multiplicity 



mg a cone size of AR < 0.2, where Ep™ included all 
ID tracks in the cone except the muon track. To help 
ensure that the muon is prompt it was required that the 
transverse impact parameter of the track do an d its un- 
certainty cr(do) satisfied \do/<r(do)\ < 3. Also the lon- 
gitudinal impact parameter Az was required to satisfy 
\Az\ < 10 mm to reduce contributions from in-time pile- 
up and cosmic ray muons. These impact parameters were 
measured with respect to the primary vertex. Events 
were rejected if there was a second muon passing the 
same kinematic selections and isolation requirements as 



above. These muon selection criteria are similar to those 
applied in Ref. [9]. 



1. Muon channel background estimates 

For the muon channel, the main backgrounds arise 
from semileptonic decays of heavy flavor hadrons in mul- 
tijet events, other leptonic decays of heavy gauge bosons, 
and tt production. The backgrounds from gauge bosons 
include W — > tv where the tau decays to a muon, 
Z — > fifi where one muon is not identified, Z — > tt, and 
diboson production. For low jet multiplicities the largest 
backgrounds are W — > tv and Z — > while for higher 
multiplicities tt production dominates (tt — > bbqq'fiv). 
Similarly to the electron channel, the number of leptonic 
background events surviving the selection criteria was es- 
timated with simulated event samples described in Sec- 
tion IV. Pythia was used only for inclusive production 
of W — > tv and Z — > tt and Alpgen for the other vector 
boson samples. The simulated leptonic background sam- 
ples were normalized to the integrated luminosity of the 
data using the predicted NNLO, NLO+NNLL or NLO 
cross sections. Discussion of the tt background follows in 
Section VE. 
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FIG. 4. Result of the £'™ lss template fits used to obtain an estimate of the multijet background for W —t fiu events with relaxed 
kinematic requirements, rriT(W) > 25 GeV and _B™ 1SS > 15 GeV. Results are shown in bins of exclusive jet multiplicity. In this 
case the multijet template was obtained with a reversed requirement on the significance of muon's impact parameter. The data 
with > 5 jets are not used for measurements because of the low event count and a poor signal-to-background ratio. 



The multijet QCD background in the muon channel 
is dominated by leptonic decays of bottom or charm 
hadrons in jets where the hadron decay involves a muon 
and neutrino. The number of background events was es- 
timated by fitting, for each exclusive jet multiplicity, the 
E™ lss distribution in the data (with relaxed selection re- 
quirements on E™ iss and m T (W): E™ iss > 15 GeV and 
m^(W) > 35 GeV) to a sum of two templates: one for 
the multijet background and another which included sig- 
nal and the leptonic backgrounds. The fit determined 
the relative normalization of the two templates. The 
shapes for the second template were obtained from sim- 
ulation and their relative normalization was fixed to the 
predicted cross sections. The full kinematic selection, 
E™ iss > 25 GeV and m T (W) > 40 GeV, was imposed 
on the multijet background samples to convert their nor- 
malization coefficients from the relaxed to full selection. 

The template for the multijet background was obtained 
from data by applying all the standard muon selection 
requirements, except that the requirement on the signifi- 
cance of the transverse impact parameter was reversed to 
\do/cr(do)\ > 3. In addition, the impact parameter was 
required to be within 0.1 < |drj| < 0.4 mm. The lower 
cut on the impact parameter reduces signal W — > fiv 



events leaking into the background sample. The up- 
per cut on | do | was placed to minimize bias from mul- 
tijet events where an isolated muon is accompanied by 
a nearby energetic jet; the isolated muons from decays 
of heavy hadrons tend to have large impact parameter. 
The background events with a muon and an energetic 
jet do not survive the standard muon selection due to 
the stringent requirement on the impact parameter, in 
conjunction with the isolation cut. 

The comparisons of the template fits to the E™ lss dis- 
tributions are presented in Fig. 4 for W — > \w events with 
the relaxed selection requirements on £"™ lss and tot(VF). 
Figure 5 shows the final mx(W) distributions in the var- 
ious bins of inclusive jet multiplicity for events passing 
the normal selection requirements. 

Another set of templates for the multijet background 
was obtained using a simulated dijet sample from 
Pythia where the event record was required to contain 
at least one muon with pt > 8 GeV. The second set of 
templates was fitted to data in the same manner as the 
first in order to estimate a systematic uncertainty in the 
number of multijet background events. The uncertainty 
increased with the jet multiplicity from 15% for the in- 
clusive VF-boson sample up to 76% for events with a W 
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FIG. 5. Comparison of transverse mass distributions rriT(W) for W —¥ [iv events. Results are shown in bins of inclusive 
jet multiplicity for events passing the normal selection requirements. MC predictions for the W — > fiv signal and leptonic 
backgrounds are normalized to luminosity using (N)NLO cross sections and the multijet background is estimated from data. 



boson and four or more jets. 



2. Muon channel systematic uncertainties 



TABLE IV. Summary of relative systematic uncertainties as- 
sociated with the muon channel. 



Quantity 


Uncertainty 


Trigger efficiency 
Muon reconstruction 
and identification 
Muon pt scale 
Muon pt resolution 
Multijet QCD background 


0.6-0.7% a 
~l.I% b 

~0.4% a 

< 6% c 

15 — 76% d ; difference between the 
two templates, see Section VD 1 



a n — pt dependent 
b n — tj> dependent 

c n — pt dependent relative to the measured resolution 
d varies with jet multiplicity 



The muon trigger efficiencies were measured using a 
Z [a/j, sample triggered by a muon candidate other 
than the muon under study [16]. Scale factors close to 
unity, relative to the value predicted by the Monte Carlo 



simulation, were obtained for the muon triggers. The 
scale factors were calculated as a function of muon r\ and 
Pt- The same sample of events was used to determine 
the muon reconstruction and identification efficiencies as 
a two-dimensional matrix in r\ and 4> [50, 51]. The mea- 
sured efficiencies were used to correct the simulated sam- 
ples. The average efficiency correction is consistent with 
unity within a systematic uncertainty of 1.1%. 

The measured momentum scale and resolution for the 
muons were studied with Z — > events [52]. The muon 
transverse momentum and its resolution were calibrated 
as a function of r\ and pt- The systematic uncertainties 
for the muon channel are summarized in Table IV. 



E. Detector-level comparisons between final states 
of W — > ev + jets and W — > \w + jets 

Observed and expected distributions for several vari- 
ables have been compared for the electron and muon 
channels. The observed distributions are shown with sta- 
tistical uncertainties. The expected distributions are pre- 
sented with experimental uncertainties that include those 
described later in Section VG in addition to the uncer- 
tainties specific to the two channels from Sections V C 2 
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FIG. 6. The uncorrected inclusive jet multiplicity distribu- 
tion. The following remarks apply to this and subsequent 
figures. Top: electron channel. Bottom: muon channel. The 
signal and leptonic backgrounds are shown using simulations, 
whereas the multijet background uses the method described 
in the text. The signal and leptonic backgrounds are normal- 
ized to the predicted cross sections. The black-hashed regions 
illustrate the experimental uncertainties on the predicted dis- 
tributions. 



FIG. 7. The uncorrected distribution in pt of the jet with 
the highest pt, in events with one or more jets. 
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FIG. 8. The uncorrected distribution in pt of the jet with FIG. 9. The uncorrected distribution in pt of the jet with 
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FIG. 12. The uncorrected distribution in y(l) — y (first jet), 
rapidity difference between the lepton and the leading jet, 
for events with one or more jets. 



FIG. 13. The uncorrected distribution in y(£) + y (first jet), 
sum of rapidities of the lepton and the leading jet, for events 
with one or more jets. 
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FIG. 14. The uncorrected distribution as a function of 
Ai?(first jet, second jet), distance between the first two jets, 
for events with two or more jets. 



FIG. 15. The uncorrected distribution as a function of 
Ac/>(first jet, second jet), azimuthal separation between the 
first two jets, for events with two or more jets. 
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first two jets, for events with two or more jets. 



and VD2. Distributions of the inclusive jet multiplicity 
are shown in Fig. 6. Figures 7, 8, 9, and 10 show dis- 
tributions in pt of the first four (highest pr) jets. The 
rapidity of the first jet is shown in Fig. 11. The difference 
and sum of the rapidities of the lepton and the first jet are 



shown in Figs. 12 and 13, respectively. Variables depen- 
dent on the azimuthal and rapidity separations between 
the first two jets are featured in Figs. 14, 15, and 16. 
Overall, a good agreement is seen between measured and 
predicted distributions. Minor discrepancies appear for 
jet pairs with large rapidity separation in Figs. 14 and 16. 
Fig. 12 illustrates discrepancies for events with the first 
jet separated in rapidity from the lepton. Predictions in 
Figs. 11, 12, and 13 are found to be sensitive to the choice 
of PDF. 

Top quark pair production is a substantial background 
to VF+jcts in events with four or more jets as can 
be seen in Fig. 6. The predicted ti cross section of 
165^g pb [24] is fully consistent with the measured 
value of 171±20(stat.)±14(syst.)l|(lum.) pb, obtained 
with the same 2010 data sample [53]. Here the predicted 
one was used to obtain the cross section results. 

Several kinematic distributions were used to check the 
normalization of the ti component in the channels with 
a W boson plus four or more jets. These included the 
rapidity of the charged lepton and the mass of the W- 
jet system. The normalizations obtained were consistent 
with the expected value but had a statistical uncertainty 
too large to usefully constrain the ti cross section. 

F. Unfolding of efficiency and resolution effects 

The yield of signal events was corrected back to the 
particle level separately for the two lepton channels, tak- 
ing into account detector acceptance and reconstruction 
efficiency. The correction was made using an iterative 
Bayesian method of unfolding [54] . Bin sizes in each his- 
togram were chosen to be a few times larger than the res- 
olution of the corresponding variable. Migration matrices 
were computed using the Alpgen W+jets event gener- 
ator plus full detector simulation, restricting the events 
to the common phase space: 

• > 20 GcV (£ = electron or muon), 

• 1^1 < 2.5, 

• p^ > 25 GeV, 

• m T {W) > 40 GeV, 

• pi^ > 30 GcV, 

• \yi ct \ < 4.4 and AR(£,jet) > 0.5. 

The common phase space requirements were applied to 
generated objects before the detector simulation. In this 
analysis, particle level jets were constructed in simulated 
events by applying the anti-fc t jet finder to all final state 
particles with a lifetime longer than 10 ps, whether pro- 
duced directly in the pp collision or from the decay of 
particles with shorter lifetimes. Neutrinos, electrons, and 
muons from decays of the massive W bosons were not 
used for the jet finding. Final state QED radiation differs 
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for electrons and muons, and its effects were corrected in 
the combined cross sections. Fiducial cross sections for 
each channel were defined using final-state leptons for 
which collinear radiation in a cone of R = 0.1 is added 
to the lepton four- momentum [55] . This accounts for the 
most significant effects of collinear QED radiation. A 
residual correction for large-angle radiation outside this 
cone is then applied to bring both electrons and muons 
to the Born level for the combined cross sections. These 
correction factors range from 0.985 to 0.995 and are sim- 
ilar for both electrons and muons. 

Instead of inverting the migration matrix, the unfolded 
distributions were determined using Bayes' theorem to 
recalculate the particle level distributions from the de- 
tector level distributions. The unfolded values were cal- 
culated using different numbers of iterations for differ- 
ent bins of a distribution. The standard Bayesian ap- 
proach treats all bins using the same number of iter- 
ations. Fewer iterations were performed for bins with 
few events than for bins with large numbers of events 
to avoid large statistical fluctuations in the tails of the 
distributions. The number of iterations was limited for 
a bin once the statistical uncertainty becomes substan- 
tially larger than the change due to the last application 
of the unfolding matrix [56] . Tests with simulated data 
showed that the iterative Bayesian method was sufficient 
to recover particle-level distributions. The dominant de- 
tector to particle level corrections in the electron channel 
come from electron reconstruction efficiency («30% cor- 
rection). In the muon channel, the dominant corrections 
come from trigger and reconstruction efficiency (correc- 
tions of ps 10 - 20% and ps 10% respectively). The sta- 
tistical uncertainty on the unfolding was estimated using 
toy simulations. The systematic uncertainties on the un- 
folding included the uncertainty on the migration matrix 
which was estimated by using the alternative Sherpa 
simulation for W+jets production (see Table I). 



G. Overall Systematic Uncertainties 

In addition to the systematic uncertainties specific to 
the electron and muon channels documented earlier in 
Sections VC2 and VD2, respectively, there are a num- 
ber of common sources of uncertainty. As a brief re- 
minder, the uncertainty on the identification efficiency for 
electrons results in ^43% variation of the Nj Ct > 1 cross 
section, giving the largest variation among the electron- 
specific uncertainties. Similarly, the uncertainty on re- 
construction and identification efficiency of muons cor- 
responds to a variation of ±1.1% in the Aj e t > 1 cross 
section and represents the single largest muon-specific 
uncertainty. 

The dominant source of systematic uncertainty in the 
cross section measurement for both electron and muon 
channels is the uncertainty in the jet energy scale [45]. 
For Ajot > 4, uncertainties on the predicted tt cross sec- 
tion and tt shape also become significant and can be as 



high as 10% and 21%, respectively. The luminosity un- 
certainty enters primarily through the signal normaliza- 
tion but also has a small effect on the estimation of the 
leptonic backgrounds. 

Uncertainties in the jet energy scale (JES) and jet en- 
ergy resolution (JER) were determined from data and 
simulation [45]. The JER uncertainty was 10% of the 
jet energy resolution [45] . The JES uncertainty varies as 
a function of jet pr and 77, and ranges from ~ 2.5% at 
60 GeV in the central region to ~ 14% below 30 GeV 
in the forward regions; the uncertainty increases mono- 
tonically with the absolute value of jet pscudorapidity. 
The uncertainty on the correction of the JES for pile-up 
pp interactions is less than 1.5% per additional interac- 
tion for jets with pt > 50 GeV. To take into account the 
differences in calorimeter response to quark- and gluon- 
initiated jets, the uncertainty on the fraction of gluon- 
initiatcd jets, the flavor composition [45], was estimated 
by comparing the fractions in Sherpa and Alpgen sim- 
ulations for W^+jets production. For jets accompanied 
by a second jet within AR < 0.7, an additional uncer- 
tainty is added to the JES uncertainty; the additional 
uncertainty is less than 2.8%. To estimate the impact of 
the JES uncertainty, jet energies in the simulated events 
were coherently shifted by the JES uncertainty, and the 
fimiss vec t or was recomputed. In addition, simulated en- 
ergy clusters in the calorimeters not associated with a 
jet or electron, such as those coming from the underly- 
ing event and pile-up interactions, were scaled using a px 
and |ry| dependent uncertainty [16], ranging from ±5.5% 
for central clusters at pt — 500 MeV to ±3% at high pt- 
Similarly the simulated jet energies were smeared by the 
JER uncertainty and the E™ lss vector was recomputed. 
The full analysis was repeated with these variations, and 
the cross sections were recomputed; the change in the 
cross section was taken as the systematic uncertainty. 
The uncertainty on the measured cross sections caused 
by the uncertainties on the JES and cluster energy scale 
increases with jet multiplicity from 9% for Nj Ct > 1 to 
37% for TVjct > 4. The impact of the JES uncertainty 
is amplified for events with high jet multiplicities due 
to the large subtraction of tt events, corresponding to 
^54% of these events. The simulated jet multiplicity of 
the top background is sensitive to the JES. The magni- 
fication is somewhat smaller when jets are selected with 
pjf > 20 GeV instead of 30 GeV; the JES-related uncer- 
tainty on the A'jet > 4 cross section is up to 29%. 

The uncertainty due to jets originating from pile-up 
interactions and the influence of the JVF selection re- 
quirement includes the efficiency of the requirement and 
how well the rate of pile-up jets is modeled in the simu- 
lation. As a conservative estimate, the percentage of jets 
in the data removed by the JVF requirement is applied 
as the uncertainty. This results in a 1.5% uncertainty for 
jets with px < 40 GeV with a resulting uncertainty on 
the cross section of 1% for Nj et > 1. 

Other uncertainties which were considered include the 
jet reconstruction efficiency and biases in the procedure 
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TABLE V. Summary of systematic uncertainties on the cross sections. The uncertainties are shown for 7Vj ct > 1 and iVjct > 4. 
The sign convention for the JES and lepton energy scale uncertainties is such that a positive change in the energy scale results 
in an increase in the jet or lepton energy observed in the data. 



W — ¥ eu channel 







Cross Section Uncertainty (%) 


Effect 


Range 


ATjet > 1 


JVjet > 4 


Jet and cluster energy scales 


2.5-14% (dependent on jet r\ and pt) 


±9.0,-6.6 


±37, -35 


Jet energy resolution 


~10% on each jet (dependent on jet ri and »t) 


±1.6 


±6 


Electron trigger 


±0.5% 


±0.6,-0.5 


±1 


Electron reconstruction 


±1.5% 


±1.7,-1.6 


±4 


Electron identification 


±2-8% (dependent on electron r\ and pt) 


±4.3,-4.0 


±10,-9 


Electron energy scale 


±0.3-1.6% (dependent on rj and pt) 


±0.6 


±1,-3 


Electron energy resolution 


< 0.6% of the energy 


±0.0 


<1 


Pile-up removal requirement 


~ 1.5% in lowest jet pt bin 


±1.1 


±3 


Multijet QCD background shape 


from template variation 


±0.7 


±11 


Unfolding 


Alpgen vs. Sherpa 


±1.5 


±6 


Luminosity 


±3.4% 


±3.8,-3.6 


±9,-8 


NNLO cross section for W/Z 


±5% 


±0.2 


<1 


NLO cross section for ti 


t!o% 


±0.3 


±10 


Simulated ti shape 


from samples with more or less ISR 


±0.1 


±12,-21 



W —J- \iv channel 







Cross Section Uncertainty (%) 


Effect 


Range 


ATjet > 1 


Ajot > 4 


Jet and cluster energy scales 


2.5-14% (dependent on jet r\ and pt) 


±8.2,-6.2 


±33, -26 


Jet energy resolution 


10% on each jet (dependent on jet r\ and pt) 


±1.5 


±5 


Muon trigger 


±0.7% (±0.6%) in barrel (endcap) 


±0.6 


±1 


Muon reconstruction 


±1.1% 


±1.1 


±2 


and identification 








Muon momentum scale 


±0.4% 


±0.2,-0.3 


<1 


Muon momentum resolution 


±6% 


±0.1 


<1 


Pile-up removal requirement 


~ 1.5% in lowest jet pt bin 


±1.0 


±3 


Multijet QCD background shape 


from template variation 


±0.8 


-20 


Unfolding 


Alpgen vs. Sherpa 


±0.2 


<1 


Luminosity 


±3.4% 


±3.7,-3.5 


±7 


NNLO cross section for W/Z 


±5% 


±0.4 


<1 


NLO cross section for ti 


±Io% 


±0.4,-0.3 


±10,-7 


Simulated ti shape 


from samples with more or less ISR 


<0.1 


±13,-15 



for correcting for detector effects (by comparing correc- 
tion factors obtained with Alpgen to those obtained 
with Sherpa). Their effect on the cross section was 
found to be smaller than the uncertainties described be- 
fore. All of these systematic uncertainties were also ap- 
plied to the estimates of the multijet and leptonic back- 
grounds in both electron and muon channels. In addi- 
tion, for the leptonic backgrounds the uncertainty in the 
NNLO cross sections was taken to be 5% for W/Z pro- 
duction as in Ref. [16]. The ti cross section uncertainty 
was taken to be 1l[ % [24]. The uncertainty on the shapes 
of the ti distributions was estimated using AcerMC sim- 
ulations where rates of ISR and FSR were altered with 
respect to the default settings. Samples with altered ISR 
were used to estimate the shape uncertainty since their 
impact on measured cross sections was the largest among 
these samples. The procedure has been used for ATLAS 
measurements involving top pair production [53]. 

The systematic uncertainties in the cross section mea- 



surement are summarized in Table V for A?j e t > 1 and 
A/j e t > 4; most of the uncertainties are approximately 
independent of the jet multiplicity, except for the uncer- 
tainty due to the jet energy scale and resolution, multijet 
background shape, ti production, and pile-up jet removal. 
The uncertainty due to the jet energy scale dominates for 
events with at least one jet as illustrated in Fig. 17. 



In the cross section ratio measurement, 
a(W+ > N ict )/a(W+ > N ict - 1), the uncertainty 
due to the jet energy scale uncertainty remains the 
dominant effect, amounting to approximately 5-20% 
on the ratio. The luminosity uncertainty does not 
completely cancel in the ratio because the background 
estimates are affected by the luminosity uncertainty 
and the background levels vary as a function of jet 
multiplicity. 
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FIG. 17. Systematic uncertainties on the cross section as a 
function of the inclusive jet multiplicity. The uncertainty due 
to the jet energy scale is bounded by the two black lines. 
The quadratic sum of the other systematic uncertainties is 
presented as the shaded area. The uncertainties are for the 
sum of the electron and muon cross sections. 



VI. NEXT-TO-LEADING-ORDER QCD 
PREDICTIONS 



The MCFM v5.8 [28] and BlackHat-Sherpa [11] 
predictions were obtained with the same jet algorithm 
and same kinematic selection requirements applied to the 
data. In both cases, renormalization and factorization 
scales were set to Ht/2, where Ht is the scalar sum 
of the j»t of all the partons and of the lepton and neu- 
trino from the JU-decay. The PDFs used for MCFM 
were CTEQ6L1 [38] and CTEQ6.6M [40] for the LO and 
NLO calculations, respectively. For BlackHat-Sherpa 
CTEQ6.6M was used for both LO and NLO calculations. 

The systematic uncertainty in the MCFM and 
BlackHat-Sherpa cross section due to renormalization 
and factorization scales were estimated by varying the 
scales by factors of two, up and down, in all combina- 
tions. The ratio of one scale to the other was kept within 
the range 0.5 to 2.0 to avoid the effects of large logarithms 
of the scale ratios in some kinematic regions. The cross 
section ratio, a(W+ > Nj et )/cr(W+ > N$ ei — 1), was re- 
calculated for each variation of the scales and the result- 
ing uncertainty was determined using the recalculated 
values. Overall, the asynchronous variations of scales re- 
sulted in bigger deviations from the nominal values than 
the synchronous variations. The upper and lower uncer- 
tainties were taken as the maximum deviations from the 
nominal value. 

Following the PDF4LHC recommendations [57], PDF 
uncertainties were computed by summing in quadrature 



the dependence on each of the 22 eigenvectors character- 
izing the CTEQ6.6 PDF set; the uncertainty in a s was 
also taken into account. The uncertainties were scaled to 
a confidence level (C.L.) of 68%. Two alternative PDF 
sets, MSTW2008 [19], with its set of 68% C.L. eigenvec- 
tors, and NNPDF2.0 [58], were also examined. The error 
envelope of CTEQ6.6 was found to contain nearly all 
variations due to the two alternative PDF sets. The un- 
certainties due to the scale variations were substantially 
larger than those due to PDFs. 

As a cross-check, cross sections from BlackHat- 
Sherpa and MCFM were compared for events with up 
to two jets, and found to be nearly identical. Therefore, 
only distributions from BlackHat-Sherpa were com- 
pared to the measured cross sections. 

Bin-by-bin corrections for non-pQCD effects, 
hadronization and underlying event, were computed 
using simulated VF+jets samples for each predicted 
distribution for the NLO cross sections. The corrections 
were taken to be the ratios of the distributions for 
particle-level jets to the distributions for parton-level 
jets, where the sample for parton-level jets was produced 
with the underlying event turned off. To calculate 
the central values, samples from Alpgen v2.13 were 
showered with Herwig v6.510 and Jimmy v4.31 set to 
the AUET2 tune [59]. The systematic uncertainty on 
the non-pQCD corrections was evaluated by comparing 
the central values to corrections from samples where 
Alpgen was showered with Pythia v6.4.21 set to the 
AMBT1 [44] event generator tune. The corrections 
and their uncertainties were applied to all the NLO 
predictions presented in the paper. 



VII. CROSS SECTION RESULTS 

The measured VF+jets cross sections were calculated 
in the limited kinematic region defined in Section V F. All 
cross sections were multiplied by the leptonic branching 
ratio, Br(W -> Iv). 

The cross sections for the W — > ev and W — > fiv chan- 
nels were calculated separately and then compared. The 
two sets of cross sections were found in good agreement 
within their uncorrelated uncertainties. The systematic 
uncertainties specific to the individual channels were con- 
sidered fully uncorrelated and the common systematic 
uncertainties fully correlated. Results for the electron 
and muon channels were combined using three passes 
of the Best Linear Unbiased Estimator (BLUE) tech- 
nique [60-62] . Three iterations were required to compute 
the upper systematic uncertainty, the central value, and 
the lower systematic uncertainty. The combination im- 
proved uncertainties and fluctuations in the tails of the 
measured distributions. 

Particle level expectations from Alpgen and Sherpa 
simulations as well as a calculation using BlackHat- 
Sherpa were compared to the measured cross sections. 
Pythia is shown only for selected distributions that are 
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FIG. 18. W+jets cross section results as a function of cor- 
rected jet multiplicity. The following remarks apply to this 
and subsequent figures unless specific comments are provided. 
The cross sections are quoted in the kinematic region de- 
scribed in Section VF. For the data, the statistical uncer- 
tainties are shown with a tick on the vertical bars, and the 
combined statistical and systematic uncertainties are shown 
with the full error bar. Also shown are predictions from Alp- 
Gen, Sherpa, Pythia and BlackHat-Sherpa, and the ra- 
tio of theoretical predictions to data (Pythia is not shown 
in the ratio). The distributions from Sherpa, Pythia, and 
Alpgen were normalized to the NNLO total W-boson pro- 
duction cross section. 



given as a function of corrected jet multiplicity. As 
Pythia features LO matrix element accuracy for events 
with up to one jet, it does not provide a good description 
of the data for jet multiplicities greater than one. The 
Alpgen, Pythia, and Sherpa predictions were normal- 
ized to the NNLO inclusive VF-boson production cross 
section. The version of BlackHat-Sherpa used here 
provides NLO predictions at parton level for VF-boson 
production with Aj ot < 4. No additional normalization 
was applied to the BlackHat-Sherpa predictions. 

The measured IF+jets cross sections and the cross sec- 
tion ratios are shown as a function of the corrected jet 
multiplicity in Figs. 18 and 19. The cross section is shown 
as a function of the px of the first jet for A/jet > 1 to 
^Vjot > 4 events separately in Fig. 20, the second jet for 
A/jet > 2 to A/jet > 4 events separately in Fig. 21, the 
third jet for Aj e t > 3 and Nj et > 4 events separately in 
Fig. 22, the fourth jet for N iet > 4 events in Fig. 23. The 
jets are ordered from the highest to lowest px- The dif- 



ai 
+ 



0.3 



aT 0.2 

CD 



Al 
+ 



0.1 



i r 

W->K> + jets 
■e-Data 2010,\/s=7TeV 
▼ ALPGEN 
A SHERPA 
" * PYTHIA 

BLACKHAT-SHERPA 



j Ldt=G 



=36 pb 
ATLAS 



, — + — 


1 






A 







>1/>0 



anti-k T jets, R=0.4 

pL et >30 GeV, |y jet |<4.4 

j L_i i 

>2/>1 >3/>2 >4/>3 

Inclusive Jet Multiplicity Ratio 



FIG. 19. W+jets cross section ratio results as a function of 
corrected jet multiplicity. 



ferential cross section as a function of i?r is shown for 
A/jet > 1 to Ajct > 4 in Fig. 24. Here H T is defined as a 
scalar sum over px of the lepton, neutrino (i5™ ISS ), and 
all jets in the event. Hf is often used to set the renor- 
malization and factorization scales in fixed-order calcula- 
tions and is therefore an interesting variable to compare 
between data and predictions. 

The measured Ht distribution for events with one 
or more jets is not well described by the BlackHat- 
Sherpa prediction. The prediction is calculated inclu- 
sively, at NLO, for events with a W boson and one or 
more jets: because of the limited order of the calcula- 
tion, matrix elements with three or more real emissions 
of final-state partons are not included in the calculation. 
In contrast, Alpgen, where LO matrix-element terms 
with up to five final-state partons are utilized, describes 
the data well. The data themselves are, as stated above, 
inclusive of all higher jet multiplicities. A modified treat- 
ment of BlackHat-Sherpa prediction was introduced, 
where higher-order NLO terms with two, three, and four 
real emissions were also added to the Aj 0t > 1 distribu- 
tion: this is shown in Fig. 25. The higher-order terms 
were combined by matching them exclusively in jet mul- 
tiplicity by counting parton jets with px > 30 GeV. The 
matching scheme is required to reduce double-counting of 
cross sections. This case illustrates the challenges of com- 
paring NLO calculations to complex inclusive jet vari- 
ables like Ht- In Fig. 26 the cross sections are shown as 
a function of the invariant mass, m,(jets), of the first two, 
three, and four jets for events with Aj et > 2, Aj Ct > 3, 
and Ajct > 4, respectively. The invariant mass of the 
multijet system is also considered for the renormaliza- 
tion and factorization scales in fixed-order pQCD calcu- 
lations. Overall, these distributions constitute a set of 
tests for factorization and rcnormalization scales used in 
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FIG. 20. W+jets cross section as a function of the pt of the 
first jet in the event. The pt of the first jet is shown separately 
for events with > 1 jet to > 4 jet. The > 2 jet, > 3 jet, and 
> 4 jet distributions have been scaled down by factors of 
f0, 100, and 1000 respectively. Shown are predictions from 
Alpgen, Sherpa, and BlackHat-Sherpa, and the ratio of 
theoretical predictions to data for > 1 jet and > 2 jet events. 



calculations of a s ; the Alpgen samples demonstrate a 
better agreement with data than Sherpa due to differ- 
ences in the scales and PDFs described in Section IV. 

Distributions dependent on rapidities of the leptons 
and the first jet are shown in Figs. 27, 28, and 29 for 
y(first jet), y(£) - y(firstjet), and y(£) + y(first jet), 
respectively. These distributions are sensitive to PDFs 
used for calculations of LO and NLO matrix elements. 
Predictions from BlackHat-Sherpa and Sherpa were 
produced with CTEQ6.6M, a NLO PDF, while Alp- 
gen used CTEQ6L1, a LO PDF. The shape of the 
distributions from Sherpa were found to be similar to 
BlackHat-Sherpa. Alpgen gave a different descrip- 
tion of the y(£) — y(first jet) distribution. The deviations 
observed between the data and BlackHat-Sherpa at 
high jet rapidities in Fig. 27 may be caused by insuffi- 
cient knowledge of the gluon PDFs at high x. 

Lastly, distances between the first two jets are ex- 
plored in Figs. 30, 31, and 32 by defining the distance as 



^ 10- 
O 



Cl 



cl 

T3 



CO 

ro 

Q 

>-, 
i 

o 

CD 



10" 
10 

1 

10" 1 
10" 2 
10" 3 
10" 4 
10" 5 

1.5 

1 

0.5 



lets 



-I I I I I r 
i Ldt=36 pb" 1 



-is 



VJ-ttv + jets 
^3-Data 2010,Vs=7TeV 

▼ ALPGEN 

A SHERPA 

BLACKHAT-SHERPA = 

ATLAS 



Jets 



X1Q-1 



~ 4 ' e ts, x i ai 



anti-k T jets, R=0.4 
r p^OGeV, |y iet |<4.4 



+ 



: W + >2 jets 



50 100 150 

Second Jet p T [GeV] 

FIG. 21. ly+jets cross section as a function of the pt of the 
second jet in the event. The pt of the second jet is shown 
separately for events with > 2 jet to > 4 jet. The > 3 jet 
and > 4 jet distributions have been scaled down by factors of 
10 and 100 respectively. Shown are predictions from Alpgen, 
Sherpa, and BlackHat-Sherpa, and the ratio of theoretical 
predictions to data for and > 2 jet events. 



Ai?(first jet, second jet), y (first jet) — ^(second jet), and 
A0(firstjet, second jet), respectively. This set of mea- 
surements offers a test of hard parton radiation at large 
angles and of matrix element to parton shower matching 
schemes. The majority of jets are modeled via the ME 
calculation for the jet pairs with large angular separation, 
when AR and Acf> are close to ir. Collinear radiation at 
small angular separation, when AR is small, is produced 
mainly via the parton shower. Overall, Alpgen and 
BlackHat-Sherpa demonstrate good agreement with 
the data while Sherpa deviates due to the differences in 
PDFs, a s , and factorization scales. 

All distributions were also produced with the selection 
requirement on pi^ reduced from 30 GeV to 20 GeV. The 
results for the softer threshold are given in Appendix A. 
The softer threshold makes the cross sections more sen- 
sitive to the non-pQCD and experimental effects, espe- 
cially for forward jets. 

All these cross sections accompanied by the non-pQCD 
and QED corrections are available in HEPDATA. 
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FIG. 22. VK+jets cross section as a function of the pt of 
the third jet in the event. The jjt of the third jet is shown 
separately for events with > 3 jet and > 4 jet. The > 4 jet 
distribution has been scaled down by a factor of 10. Shown are 
predictions from Alpgen, Sherpa, and BlackHat-Sherpa. 
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VIII. CONCLUSIONS 

This paper presents a measurement of the H^+jets 
cross section as a function of jet multiplicity in pp colli- 
sions at y/s = 7 TeV in both electron and muon decay 
modes of the W boson, based on an integrated luminos- 
ity of 36 pb -1 . The ratios of cross sections a(W+ > 
Nj et )/(r(W+ > iVj e t — 1) have been calculated for inclu- 




200 



400 



600 

H T [GeV] 



FIG. 24. W+jets cross section as a function of Hi, shown 
separately for > 1 jets to > 4 jets. The > 2 jet, > 3 jet, 
and > 4 jet distributions have been scaled down by factors of 
10, 100, and 1000 respectively. Shown are predictions from 
Alpgen, Sherpa, and BlackHat-Sherpa, and the ratio of 
theoretical predictions to data for > 1 jet and > 2 jet events. 
The apparent discrepancy between the data and BlackHat- 
Sherpa predictions is discussed in the text. 



sive jet multiplicities, iVjet, that range between 1 — 4 for 
the > 30 GeV jet threshold and between 1 — 5 for 
the > 20 GeV threshold. Measurements are also pre- 
sented of the px distribution of the first through fourth 
jets in the event, of the invariant masses of two or more 
jets, of the distances between the lepton and the first 
jet, of the distances between the first two jets, and of 
the Ht distribution. The results have been corrected for 
all detector effects and are quoted in an ATLAS-spccific 
range of jet and lepton kinematics. This range is al- 
most fully covered by the detector acceptance, so as to 
avoid model-dependent extrapolations and to facilitate 
the comparison with theoretical predictions. Good agree- 
ment is observed between the predictions from the multi- 
parton matrix element generator Alpgen and the mea- 
sured distributions. At the same time, Sherpa demon- 
strates a slightly worse agreement with the experimental 
results than Alpgen. The paper features the first com- 
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FIG. 25. VF+jets cross section as a function of Ht, shown 
for final states with > 1 jets. The cross sections are quoted 
in the kinematic region described in Section VF. Shown are 
predictions from Alpgen, BlackHat-Sherpa, and modified 
BlackHat-Sherpa and the ratio of these NLO theoretical 
predictions to Alpgen. The BlackHat-Sherpa predictions 
were modified by introducing higher-order NLO terms with 
two, three, and four real emissions to the 7Vj ot > 1 distribu- 
tion. The distribution from Alpgen was normalized to the 
NNLO total VF-boson production cross section. 



parison between the NLO predictions and the LHC data 
for events with a W boson and four jets. Calculations 
based on NLO matrix elements in MCFM (available for 
jet multiplicities iVj et < 2) and in BlackHat-Sherpa 
(available for jet multiplicities N- Jc t < 4) are generally in 
good agreement with the data; deviations are observed 
in the dcr(W4- > jet)/df/r distribution at large Hi and 
in the tails of der/dy(jet) and dcr/d(y(^) — y(jet)) distri- 
butions. 
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FIG. 31. VK+jets cross section as a function of j/(first jet) — 
^(second jet) for events with > 2 jets. Shown are predictions 
from Alpgen, Sherpa, and BlackHat-Sherpa, and the ra- 
tio of theoretical predictions to data. 
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Appendix A: Results for a jet threshold of pt > 20 
GeV 



Here we present results for jets selected with a 20 GeV 
threshold in px- The distributions are the same vari- 
ables as for jets with the 30 GeV threshold shown in 
Section VII except that the data with iVj et > 5 was 
used for physics conclusions; the 20 GeV threshold im- 
proved the signal-to-background ratio and event count. 
The softer threshold makes the cross sections more sen- 
sitive to the non-pQCD and experimental effects such 
as the underlying event model, multiple parton interac- 
tions, parton fragmentation, hadronization, and pile-up 
pp interactions. The corrections accounting for the non- 
pQCD effects, that were applied to BlackHat-Sherpa 
calculations, increased monotonically with the absolute 
value of jet rapidity from ~1.0 up to ~2.4. The uncer- 
tainties on the corrections are also larger in the forward 
region. 
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FIG. 34. W+jets cross section ratio results as a function 
corrected jet multiplicity. 



FIG. 33. W+jets cross section results as a function of cor- 
rected jet multiplicity. The following remarks apply to this 
and subsequent figures unless specific comments are given. 
The cross sections are quoted in the restricted kinematic re- 
gion described in Section VF, except here p^T > 20 GeV. For 
the data, the statistical uncertainties are shown with a tick 
on the vertical bars, and the combined statistical and sys- 
tematic uncertainties are shown with the full error bar. Also 
shown are predictions from Alpgen, Sherpa, Pythia and 
BlackHat-Sherpa, and the ratio of theoretical predictions 
to data (Pythia is not shown in the ratio). The distributions 
from Sherpa, Pythia and Alpgen were normalized to the 
NNbO total IT-boson production cross section. 
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FIG. 35. PF+jets cross section as a function of the pt of the 
first jet in the event. The pt of the first jet is shown separately 
for events with > 1 jet to > 4 jet. The > 2 jet, > 3 jet, and 
> 4 jet distributions have been scaled down by factors of 
10, 100, and 1000 respectively. Shown are predictions from 
Alpgen, Sherpa, and BlackHat-Sherpa, and the ratio of 
theoretical predictions to data for > 1 jet and > 2 jet events. 
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FIG. 41. ly+jets cross section as a function of y (first jet) for 
events with > 1 jets. Shown are predictions from Alpgen, 
Sherpa, and BlackHat-Sherpa, and the ratio of theoretical 
predictions to data. 
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FIG. 42. W+jets cross section as a function of y(i) — 
j/(firstjet) for events with > 1 jets. Shown are predictions 
from Alpgen, Sherpa, and BlackHat-Sherpa, and the ra- 
tio of theoretical predictions to data. 



FIG. 43. W+jets cross section as a function of y(£) + 
j/(firstjet) for events with > 1 jets. Shown are predictions 
from Alpgen, Sherpa, and BlackHat-Sherpa, and the ra- 
tio of theoretical predictions to data. 
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FIG. 44. W+jets cross section as a function of 
Ai?(first jet, second jet) for events with > 2 jets. Shown are 
predictions from Alpgen, Sherpa, and BlackHat-Sherpa, 
and the ratio of theoretical predictions to data. 



FIG. 45. VK+jets cross section as a function of j/(first jet) — 
^(second jet) for events with > 2 jets. Shown are predictions 
from Alpgen, Sherpa, and BlackHat-Sherpa, and the ra- 
tio of theoretical predictions to data. 



33 



o. 



g 
To 

T3 



CO 

ro 
Q 

>-, 
k_ 
o 

CD 



300 

250 

200 

150 

100 

50 
1.5 



-. i i i | i i 
W->k> + jets 
I -e- Data 201 0,\/s=7TeV 
▼ ALPGEN 
A SHERPA 

BLACKHAT-SHERPA 



anti-k T jets, R=0.4 
; p i r et >20GeV, |y iet |<4.4 



0.5, 



I 



Ldt=36 pb 



ATLAS 



W + >2 jet 



1* 



— T 



A A A 



A A A A 



1 2 3 

A(j)(First Jet,Second Jet) 
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Ac/>(first jet, second jet) for events with > 2 jets. Shown are 
predictions from Alpgen, Sherpa, and BlackHat-Sherpa, 
and the ratio of theoretical predictions to data. 
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